The dynamical role of initial correlation in the exactly solvable dephasing model 
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We investigate the effects of the initial correlation on the dynamics of open system in the ex- 
actly solvable pure dephasing model. We show that the role of the initial correlation come into 
play through a phase function and a weight factor, which would perform oscillations during time 
evolution, and find that the decoherence of a qubit coupled to a boson bath is more enhanced with 
respect to a spin bath in the short time. We also demonstrate that the trace distance between two 
states of a qubit can increase above its initial value, and that the initial correlation can provide 
another resource for the damply oscillation and revival of the entanglement of two qubits. We fi- 
nally investigate the dependence of the crossover of decoherence from the dynamical enhancement 
to suppression under the bang-bang pulse control on the initial correlation and the statistics of the 
bath constituents. 



PACS numbers: 



I. INTRODUCTION 



The conventional method of studying open quantum 
system is the quantum master equation in the Lindblad 
form, which relies on the uncorrelated initial condition 
and the following two assumptions The first is that 
the coupling of the system-bath interaction is weak, and 
the other is that the relaxation time of bath is much 
shorter than the response time of system. Then the sec- 
ond order perturbation is applied and the memory effect 
of the bath can be neglected for simplification. How- 
ever, these assumptions can be violated, where the pro- 
jection operator technique can provide the generalized 
exact master equation for treatment Usually it is 
impossible to find a general scheme to solve for most of 
open quantum systems. 

Even the open system initially cannot be completely 
isolated from bath, the uncorrelated initial condition is 
still useful in many cases when the coupling between the 
system and bath is weak and when we are only inter- 
ested in the long time behavior of system. In such case, 
the evolution of the system can be simply described by 
the completely positive map @ , which leads to many pro- 
found results. For example, the trace distance between 
two states of the open system can not increase above its 
initial value [!, H| . 

On the other hand, the assumption of the initially un- 
correlated condition is too restrictive and unrealistic. An 
alternative to open quantum systems is based on the ex- 
actly solvable model, which has been used to construct 
the corresponding exact master equation [5| andpinpoint 
the limitation of second order approximation Q. Based 
on some solvable models, it has been discussed that the 
initial condition can significantly influence the time de- 
velopment of system. In the pure dephasing model with 
the correlated initial condition, the decoherence function 
can have some sharp peaks in comparison with the un- 
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correlated initial condition Q. In Ref. Q it was shown 
that the exact master equation of the independent oscil- 
lator model with the uncorrelated initial condition can 
lead to some divergence problems and the negativity of 
the density matrix, which can be consistently avoided by 
the correlated initial condition. 



In this paper we address the influence of the initial 
correlation on the dynamics of the open system in the 
exactly solvable pure dephasing model. We show that 
the role of the initial correlation comes into play through 
a phase function and a weight factor, which take into 
account the memory effect of bath and might give rise 
to coherence oscillation during the time evolution. In 
particular, for the short time, the decoherence of a qubit 
coupled to a boson bath is more enhanced with respect to 
a spin bath. We also demonstrate that the trace distance 
between two states of a qubit can increase above its initial 
value, and that the entanglement of two qubits, each lo- 
cally interacting with an independent bath, can damply 
oscillate and revive to a large amount being compara- 
ble to the uncorrelated initial condition [8J. Moreover, it 
is known that the decoherence can be effectively inhib- 
ited with the bang-bang control pulses Hand depends 
on the statistics of the bath constituents Q. We finally 
investigate how the initial correlations affect the decoher- 
ence and its crossover from the dynamical enhancement 
to suppression under control. 



The organization of this paper is as follows. In Sec. II 
we briefly review the main features of two pure dephas- 
ing models with boson and fermion baths respectively. 
The dynamics of the two models without initial correla- 
tions are analyzed in Sec. Ill, including the bang-bang 
control pulses. We then study the influence of the ini- 
tial correlation on the dynamics of the qubit in Sec. IV. 
The numerical results are given in Sec. V. At last, our 
conclusion is given in Sec. VI. 



II. MODELS FOR PURE DEPHASING 



can be simplified as 



The Hamiltonian describing a two-state system (S), 
i.e. qubit coupled to a bath (B) for pure qubit-dephasing 
in the boson-bath model can be written as 0, [l(| 



H 



Ha + H f 



= u Q S z + ^uj k ala k + S z ^g k (al + a k ), (1) 



where ujq is the excited energy of the qubit, and S z is the 
spin matrix of z component with basis 2S' Z |±) = ±|±). 
The annihilation and creation operators a k and a k are 
the bath mode with frequency u k . The coupling strength 
between the qubit and the fcth bath mode is denoted by 
g k . Suppose the initial state of the system at t — is 
given by a density matrix p(0), then the state at time t 
is 



P (t) = e- iHt P (oy 



(2) 



The reduced density matrix of the qubit is the partial 
trace taken over the bath modes of the total density ma- 
trix, p s = Tib p. 

It would be convenient to work in the interaction pic- 
ture, where the new state and time evolution are given 
by 



J{Hs+H B )t 



Q(t) 

U(t,to) = Texp 
H(t) = S z ^g k ( 



p(t)e 



-i(Hs+H B )t 



i I dsH(s)}, 

to 

e iuJkt at + e 



-iuikt 



a k ). (3) 



Here the notation T represents the time-ordered product. 
The evolution of the new density matrix can be written 
as 



g(t) =U{t,t Q ) g (t )U i (t,t Q ) 



(4) 



In order to further simplify U(t,ta), we use the general- 
ized Baker-Hausdorf formula [Tj, 



Texp 



x exp 



i l dsA(s) 

to 

I '* 
2 



exp 



dsA(s) 



d Sl / ds 2 [A( Sl ),A(s 2 )] } , (5) 

to J to 



which is satisfied if the commutator [A(si), A(s2)} is a 
c— number. For the interaction Hamiltonian T-L(s) we 
have 



[H(t),n(t')} = - l -Y,glsmu k {t-t') 



(6) 



With Eqs. ([3]) and ([5]), the evolution operator U(t,to) 



U(t,to) = U + (t } t )\+)(+\ + U-(t,t )\-)(-\, (7) 
where the unitary operators U± are 



-C k {t-h)e-^>a k ] 



(8) 



The functions f(t) and £, k (t) are determined by the cou- 
pling to the bath, 



= ^f^Kt-sinwfei), 



fc 

&(t) = ^(l-e*"*). 
Wfc 



(9) 



Similarly, the Hamiltonian for a qubit coupled to a 
spin-bath is given by @ 

H = Hs + H B + H int 

= u) S z + J2"k^- + S z J29k(<r k +v k ), (10) 
fc fc 

where the coupling strength g k is identical to that of Eq. 
([T]), and u> k denotes the excited energy of the fcth bath 
spin <7fc. It is worth to point out that in Q a factor of 1/2 
is missed in Hb of the above Hamiltonian, which made 
the authors of [6| draw incorrect conclusions that a spin 
bath leads to more harmful effects than a boson bath and 
the qubit dephasing is more strongly suppressed in a spin 
bath than in a boson bath as the bang-bang control pulse 
interval becomes shorter. This spin-bath model can be 
solved by expressing the Hamiltonian H as 



H = H +\+)(+\+H-\-)(- 

ujo 1 
2 2 



= ±^ + i$> fe 4±^). (11) 



The unitary evolution operator U(t) = e is given by 

U(t) = U + \+)(+\ + U-\-)(-\, 
u± = e =F<wot/aj-J^± i ( 12 ) 
fc 

where the unitary U k is 

Qkt , ( ' Wfc z g k x \ . $l k t 



(13) 



U? = cos— -i^-i^Jsin^- 



^fc = v^i + ffl 



It can seen from Eqs. (|7|) and (|T2|) that the popula- 
tions of the |±) states, p~s + (t) and pg~ (t) do not change 
with time, namely the pure dephasing model can not de- 



scribe the relaxation process of the qubit to equilibrium. 
However, the exact solutions of these simple models with 
initial correlations still capture the essential features of 
decoherence in the more complicated dissipative models. 



III. DEPHASING WITHOUT INITIAL 
CORRELATIONS 



The uncorrelated initial condition means that the ini- 
tial state of the total system is a direct product of the 
qubit state and the bath state, 



p(0) = ps(0)®p B 
ps{t) = Tt b [Ups(0)®pbU^]. 



(14) 



Note the map from ps(0) to ps{t) is completely posi- 
tive under the uncorrelated initial condition. Usually, 
the initial bath is assumed to be in the thermal equilib- 
rium at temperature T = 1//3, Pb — eT^ HB JZb with 
the normalization Zb = Trs [e~ ,3H£! ]. Because the evo- 
lution is governed by the unitary operator of the form 
C/+|+)(+| + U-\— )(— |, the coherence of the qubit pg~ 
evolves in time as 



P + S -(t) 

r«(t) 



-In (Ti b [U +Pb u!_] 
-In (tt b [pb UlU+] 



(15) 



where T a (t) denotes the decoherence function for the 
boson-bath model. Using Eq. ([8]) and the identity 

Tr B [p B exp ^2 ^4 ~ Vtai^j } 



=«p -E^h oth T)< (16) 



we obtain 



Ulu + = exp{£[&(*)4-£(*K] 

^ k 

r «W = 2^5fcC0th— 2 . (17) 



Eqs. pT|) and (fT8|) are the exact results for the decoher- 
ence functions under the uncorrelated initial conditions 
in the pure dephasing models. Unlike the boson-bath 
model, the decoherence function for the spin-bath model 
is independent of temperature and can lead to complete 
dephasing when 1 — 2(p|/0|) sin 2 (Slfc£/2) = at certain 
time if one of the coupling strength g k > Wfe. This is 
quite different from Eq. (|17[) in the boson-bath model, 
which does not have complete dephasing even if all the 
couplings satisfy g k > uj k . Moreover, due to the renor- 
malized spin-bath mode frequency Q k > Wk> it induces 
faster dynamics than boson-bath at zero temperature de- 
spite the same bath coupling spetrum. 

In the short time case t — > 0, we see that both Eqs. 
(|17[) and (|18p perform the quadratic behavior, namely 
T(t) cx t 2 . On the other hand, in the weak coupling limit 
cjk <C LOk and flk ~ wjt, we can see Eq. (fT5]) approaches 
to (fT?) at zero temperature since — ln(l — x) ~ x for 
x <C 1. It is known that for the boson-bath model the 
non-Markovian master equation up to the second order 
perturbation can give the exact result for the decoher- 
ence function even in the strong coupling limit because 
the expansion of the evolution unitary operator can be 
truncated at the second term, which is not legitimate for 
the spin-bath model. The quantum master equation in 
the second-order approximation is [l], [ll| 



dgsjt) 
dt 



ds Tr B [H(t), [H{s), 0s (t) ® g B ]]. (19) 



For the pure dephasing model, it leads to 

where the kernel for the boson bath is given by 
K.(t, s) — 9k cos[cjfc(i — s)] coth — 

k 



(20) 



(21) 



The solution of Eq. (|2"0"|) gives the decoherence function 
T<f E {t) of the form (|17p. which is valid as long as the 
second order analysis is legitimate, i.e. gk/uJk < 1. On 
the other hand, in the spin-bath model the kernel is 



k 



gl cos[uj k {t - s)]. 



(22) 



For the spin-bath model under the uncorrelated initial 
condition, we use Eqs. (|T2"j) . (IT3|) . and (|T4"|) to get 



pr(t) = 



T a (t) = -InlTrBlpslliU^U, 



The solution of Eq. J20) is then 



k k 



(23) 



which is also independent of temperature and approaches 
to (fT5|) only when gk/oJk *C 1. Therefore, the condition 
of the validity of the second order master equation for the 
spin-bath is more stringent than that of the boson-bath. 

Next, we analyze the effect of the dynamical decou- 



pling pulses along x direction at intervals r under the 
uncorrelated initial condition 0, Q . The simplest bang- 
bang pulses that can significantly reduce the dephasing 
rate is realized by frequent tt pulses along x direction ap- 
plied on the qubit. That is after such a pulse, the qubit 
states change as |±) — > |=f) in the Schrddinger picture. In 
the presence of the decoupling pulses, at time t = 2Nt, 
the evolution unitary operator becomes 

U(t,0)=U+\+){+\+U-\-){-\, (24) 

where 

U ± =e ± ^ 2 U T (t,t~T)---U T (2T,T)U ± (T,0), (25) 

and the phase factors e ±lLJot / 2 appear due to the state 
transformation from the Schodinger picture to the inter- 
action picture. Substituting Eqs. (j8]) into (|24|) and ne- 
glecting state-independent global phase factors that are 
irrelevant to the density matrix, we can treat the factors 
in Eq. (|24[) as commuting operators and get 



W± = e ±-at/ 2exp | ± l^ [%(r)a t 



rf k {r)a k ] , (26) 



where 



N 



%(t)=6(t)(1 



tT )E< 



,2i(n— 



(27) 



Repeat the steps to derive Eq. (fT7) . we have 



we have 

{U k -U+) N = (\t) N K)(v+\ + (X k ) N \v k )(v k \ (33) 

where the corresponding eigenvalues and eigenvectors are 
given by 



x k ±i\/l- x\ 



ia±\ 



a± 



Vk ± a/T 



vl 



(34) 



The eigenvalues of U k U k are the same as U k U k , and the 
corresponding eigenvectors are given by the replacements 
a± — » — a±. The final result we obtained is 



Kit) 



k 



Ps(°)> 
In (1 - 8F 2 ) , 



(35) 



F k = 



VL k T 



s\nN(j)k \ gk^k . 2 

— 5— sin — 

sin0 fc J Uf, 2 



where (f> k — cos -1 x k . Under the periodic pulse control, 
the expression for Y J (t) is still temperature- independent 
and approaches to Y^(t) in (|28[1 at zero temperature as in 
the weak coupling limit gk/^k <C 1 and cp k ~ ^fcT, which 
can also be obtained from the second order analysis. 



- P -rj (*)„+- 



Pj-(O), (28) 

E2 n t 2 UJ k T 1 - COS Ul k t 
5fcC0th _tan . 



In the spin-bath model, the evolution under periodic 
7r pulses is given by the unitary 

U{t) = (U_U + ) N \+)(+\ + (U + U.) N \-){-\. (29) 

Hence it is sufficient to find a closed form for the matrix 
{U-U+) N . Since 



(U-U + ) N = H{u-uf_ 



N 



(30) 



IV. DEPHASING WITH INITIAL 
CORRELATIONS 

One natural way to implement the initial correla- 
tion between the qubit and bath is to use the positive 
operator- value measurement (POVM) Q, from which we 
choose some operators E m acting on the whole system. 
If at first we take the whole system in the thermal equi- 
librium state p — e~P H /Tr[e _/3ff ], then after the action 
of POVM, the whole system becomes to be 



(36) 



and 



x k I ■ 



•''/, = I - 2 ^f si( > 



VkO~ k + \/ I 

2 &kT 



4 



(31) 



Uk 



k 



Uk 

— siniifc-r, 



the expression of [U k U k ] can be evaluated by writing 
U k U k in terms of its eigenvectors. That is if 



(32) 



where the factor Z is the normalization of p(0), 

Z = ^n ± Tr B [e-^±], (37) 

k k 

where the notation u± = (±\E} n E m \±) : and the fact that 
the total Hamiltonian is diagonal in the basis of S z has 
been used. 

In such way we prepare an initial state by measure- 



ment, instead of the supposed uncorrelated state Eq. 
(fT4| . In contrast tothe form of (|T4|) . the initial density 
matrix (j36j) is expressed in terms of the total Hamilto- 
nian H and takes into account the initial qubit-bath cor- 
relation through the interaction term H lrA of H. Con- 
sequently, the bath is no longer in thermal equilibrium 
initially, and its initial state becomes 

m 



Eq. (|38l) is quite different from the thermal state p B 



-PH. 



B jZ B which does not contain the interaction terms. 

Now we calculate the evolution of coherence p$~ (t) 
with the initial correlation introduced by Eq. (1551) . Using 
Eq. ((4]), we first write 

g+-(t) = Tr[a-Q(t)] 

= \lZ ^B^s[v-U{t)E m e- im EilA\t)] 

m 

= I^(±| J Bt i <7_ J B m |±)Tr fl [[/it/ + e-^±] 



^]>>±Tr B [[/l[/ +e -^±], 



(39) 



where the factor w± is w± = ^ m (±|£ , t i cr_i? m |±). As 
TT: B [UlU + e-' 3HB ] appeared in Eq. |J5J|, we need to find 
the expression for Tr s [C/iC7 + e- /3ff ±]. Note the identity 



1 



Tr B [e- fSH± e Vka ^ ak } 



-•7*0*1 



x exp 



^4"4 Wfc + 2 



(40) 



where a unitary transformation that diagnoses H ± was 
performed inside the trace function, and the phase is 



e 



k 



(41) 



Using Eqs. (fT6| and (|T7| . we get the final result for the 
coherence 



pt-(t)=pt-(0)e- rM 



J2 ± W±eT/3"o/2 

= p+- (0)e- r » W [cos 9 a (t) + i W sin 9 a (t)} , (42) 

where the decoherence function T a (t) has the same form 
of (fl7|). The phase function is given by 



e.(*) = £^ 



sin cjfet 



(43) 



and the weight factor is 



IT' 



w+e 



-/3w /2 



w+e -/9^o/2 _|_ w _ e P^o/2 ' 



(44) 



which take into account the memory effect of bath and 
the role of initial correlation. Here the initial coherence 
is given by 



/4~(o) 



w+e 



+ w_e 



-/3tL?o/2 i w _g< s "o/2 



it+e' 



(45) 



Similar results have been obtained in [7]. 
Similarly, for the spin-bath model, we get 

Ps~{t) = Ps~ (0)e- r " (t) [cos 9 CT (t) + «W sin 6 (7 (t)] 

xII^ + ^W' ( 46 ) 

where 



= J^tan ^(t), 



%(*) = 



^-tanh^sinOfet 
I-2l|sin 2 ^y 



(47) 



The memory effect of bath is taken into account by the 
temperature dependent function T)k(t), which also makes 
the decoherence function temperature independent. For 
convenience, we express Eqs. (|42l) and (|4"o) in a compact 
form, 



Ps~(t) = p^(0)e^ r(t) [cos6(i) + iWsm&(t)} 



ni>+w> 



(48) 



where the exponent s = for the boson bath and s = 1 
for the spin bath. 

Next, we consider the quantification of the degree of 
initial correlation Q. We note that the trace distance of 
any two states p\ and pi defined by D(pi,p2) = \p\ — 
P2I/2 with \A\ — TrV A^A cannot increase above its ini- 
tial value under the completely positive map, whereas the 
necessary and sufficient condition for the general dynami- 
cal map from p s (0) = Tr B [p(0)\ to p s (t) = Tr B [Up(0)W] 
being completely positive is iff p(0) is purely classically 
correlated [12j ■ So a witness for initial system- bath corre- 
lation is proposed in [3j , namely the amount of trace dis- 
tance of two states increased above its initial value. For 
example, we consider the quantity D(ps{t), ps(t)), where 
ps(t) is constructed from the fully uncorrelated marginal 
state of p(0), i.e. p s (t) = Tt b [Up(0)U^] with p(0) = 
ps(0) ® Pb(0). We have obviously D(p s (0), p s (0)) = 0, 
so it can be used as a witness for the initial correlation as 
argued in Q and an upper bounded was also found to be 
D(p(0),p(0)) > D(ps(t),p s (t)). By direct calculations, 



we obtain 



D(p s (t),ps(t)) 



\p+-{0)(W -W u )sm&\e- r 



where 



W u = 



-Puio/2 



2\ s / 2 

Vk) ■ 



/3^o/2 



_ e -/9w /2 _|_ u _ e f3u /2 ' 



(49) 



(50) 



If W — W u 7^ 0, the trace distance can take nonzero value 
as time going on, which indicates the initial non-classical 
correlation between the qubit and bath. 

We also note that the initial correlation can provide a 
new resource for the entanglement oscillation and revival. 
Suppose two qubits each locally interacting with an in- 
dependent bath and the initial state is prepared by the 
operator E m = |*)(*| with |¥) = (| + +} + | - -))/v / 2, 
which gives an "X" structure density matrix maintained 
during the evolution (8|. The entanglement measure of 
the two qubits, i.e. the Wootters concurrence, can then 
be obtained by 



C(t) = max [0,2|p 14 (t)| - 2^ p 22 (t)p 33 (t)\ , (51) 

where the standard product basis {|1) = | + +), |2) = 

I + -), 1 3) = I - +), 1 4) = I )} were used and the 

nonzero element pu{0) = 1/2. For such initial state, we 
have p 22 (t) = p3s(t) = 0, and 



Pu(t) = pu(0)e 

xB 1 



-2F 



(cos 29 + iW e sin 29) 



vl) 



(52) 



where the weight factor is W e = — tanh/3ojo- The re- 
sult for the uncorrelated initial condition is obtained by 
putting 9 = r/k = in the above equations. If we replace 
the initial state by |$) = (| + -) + | - +))/V2, the ef- 
fect of initial correlation would cancel out in the phase 
function, namely 9 = in the entanglement measure. 

The dynamical decoupling effect of the periodic pulses 
can be derived as in the previous section. The final re- 
sults have the same forms of (|42)) and (|46|) with the re- 
placements r — > T x , 9 a —> 0£, and r] k —> rj k , where 



E 2 CJfcT 1 — COS UJkt 
9k tan — 

k 



2^1 tanh ^ tan %^ sin 2 N(f) k 



1 



^ tan 2 ^f- cos 2N(j> k 



(53) 
(54) 



Obviously, in the limit of the pulse interval t — > 0, we 
have — > 0, OJ — > and — > 0, which means that 
the decoherence are inhibited significantly. For the finite 
pulse interval, the behaviors of the two models are quite 
different as discussed below. 
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FIG. 1: Time evolution of the coherence for the uncorrelated 
(thin lines) and correlated (thick lines) initial conditions with 
the Ohmic bath spectral density. The coherence in the weak 
(strong) boson & spin bath is represented by the dot-dashed 
(dashed) & dotted (solid) lines. Here the initial coherence is 
rescaled to unity. 




FIG. 2: (a) The evolution of the trace distance of two states 
of qubit. The dot-dashed (dashed) & dotted (solid) line is for 
the weak (strong) boson & spin bath, (b) The evolution of 
the concurrence of two qubits each locally interacting with an 
independent bath. 



V. NUMERICAL RESULT 

Because the decoherence function for the boson bath 
linearly depends on the coupling strength g\ , we can treat 
the bath within the classical noise description @,[T1|. By 
introducing the spectral density of the boson-bath, 

k 

/>oo 

-> g 2 Y^5(u-u k ) =g 2 / duN(u), (55) 
k Jo 

where the second line holds for the all equal couplings 
with the boson bath, and the number distribution of bath 
mode is denoted by N(w). Then the summation over the 
bath modes in the decoherence functions can be trans- 
formed into integral. For example, Eq. (I17|) becomes 



T a (t) = J dwJ(w)coth 



/3w 1 — cos uit 



cT / dioN(u}) coth 



/3u) 1 — cos wt 



(56) 



For the spin bath, it can not be written as an integral as 
opposed to the boson bath with the chosen bath spectrum 



since the mode frequency is renormalized to be ilk — 
VwJ+sf) which nonlinearly depends on g 2 except in 
the weak coupling regime [ijj. In order to compare the 
two models with closed forms, by inspecting Eq. (|55[) we 
suppose the qubit is equally coupled with each mode of 
the bath, i.e. gu = g for all k, but the number of each 
mode is distributed according to the chosen spectrum 
N(lu). For example, Eq. (TT5|) can thus be written as 



(a) 



dujN{ui) In 1 



(57) 



where the renormalized frequency is f2 = \J ui 2 + g 2 . 

For the numerical analysis, we choose the operator E„ 
to be 



E m = I + ^{<T X + <T Z ), 

and the constants u± and w± are thus given by 



U±= 2 ±h 



1 , 1 
w+ = — ± — . 

± 2 4 



(58) 



(59) 



Suppose the coupling constant g = 0.02, the tempera- 
ture j3 = 2, and the energy cjq = 0.1 [IS ]. We also set the 
number distribution of bath mode be the Ohmic spec- 
trum with a high frequency cut-off, N(ui) = AiVou;e _w ' A 
with the cut-off A = 5 and N = 2.5 x 10 3 , where the 
weak & strong bath corresponds to A = 1 & 10. 

In Fig. 1 we show the free evolution of the coherence 
in the bath with and without initial correlation. It can 
be seen that the initial correlation can lead to damped 
coherence oscillation for the strong bath, and reduce the 
coherence to some extent for the weak bath. We also 
see that for the short time, the boson bath dephases the 
qubit to more degree than the spin bath with the same 
initial condition. 

The evolution of the trace distance D(ps(t), ps(t)) is 
plotted in Fig. 2 (a), which can take nonzero value at 
t > witnessing the initial quantum correlation. Fig. 2 
(b) shows that the entanglement measure, i.e. concur- 
rence C, can revive to a large amount in contrast to the 
uncorrelated initial condition for the strong bath, which 
only appears as the two qubits are cou pled to a common 
strong bath without initial correlation [l5[ . We also note 
that both quantities can perform damped oscillations for 
the strong baths. 

Fig. 3 (a) shows the dependence of the coherence at 
time t = 2 in the strong baths on the pulses applied 
at a rate r = 1/N, from which we see that the ini- 
tial correlation always reduces the efficiency of the con- 
trol pulses relative to the uncorrelated case. Fig. 3 (b) 
shows the evolution of the relative change of the coher- 
ence S = \Ps~^/ Ps~ | under pulse control applied at a rate 
t = 0.15 with respect to its free evolution in the weak 
baths. It indicates that the n pulse control suppresses 
the decoherence for the boson bath, while the same con- 
trol can only suppress the decoherence for the spin bath 
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FIG. 3: Time evolution of the coherence in the strong baths 
for the uncorrelated and correlated initial conditions under 
the bang-bang control pulses, (a) The dependence of the co- 
herence at time t = 2 in the strong baths on the pulses applied 
at a rate r = 1/N. (b) Time evolution of the relative change 
of the coherence 5 = \Ps^/Ps~\ under pulse control applied 
at a rate r = 0.15 with respect to its free evolution in the 
weak baths. 
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FIG. 4: The crossover of the qubit coherence from the sup- 
pression to the enhancement under the bang-bang control 
pulses is plotted as a function of the pulse interval (r) versus 
the pulse number (N). The solid (dot-dashed) line is for the 
weak boson bath with(out) initial correlation. The dashed 
(dotted) line is for the weak spin bath with(out) initial corre- 
lation, where the dotted line also represents the crossover in 
the zero temperature bath. For the strong bath, the crossover 
merges into the case of the weak bath without initial correla- 
tion. 



in the short time duration, and then crossovers to the 
enhancement of the decoherence as time becomes longer. 
Moreover, the decoherence is more strongly suppressed 
in the boson bath than in the spin bath. 

In Fig. 4, we plot the crossover of the qubit coherence 
from the suppression to the enhancement under the bang- 
bang control pulses as a function of the pulse interval (r) 
versus the pulse number (N) . It is shown that the initial 
correlation affects the crossover of decoherence to some 
extent only at the beginning of evolution, and induces 
relatively more changes for the boson (weak) bath than 
for the spin (strong) bath. The underlying reason is that 
the bath's memory of the initial correlation would be 
lost as time being longer. On the other hand, the higher 
temperature for the boson bath, the more efficient for the 
control pulses to make the crossover occur. As shown in 



Fig. 4, the crossovers for the correlated and uncorrelated 
initial conditions in the limit of t —¥ oo are essentially 
coincident with each other. Asymptotically for the boson 
bath, note that 



suits shown in Fig. 4. For the spin bath, the crossover un- 
der control shows the temperature insensitivity instead. 
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where np(uj) = g N(ui) coth(/3w/2) and the following 



identities have been used [lj 
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Here the pulse interval t is fixed and t — 2Nt while the 
pulse number N — > oo. The crossover at distant time 
between the two regimes takes place at r = r* where t* 
is determined by the equation, 



4 / 7T \ 7T 

(^J = 2^ (0) ' 



(62) 



the solution of which is given by r* = at zero temper- 
ature and r* s» 0.289 at /3 = 2 same as the numerical re- 



VI. CONCLUSION 

We addressed the effects of initial correlation on the 
dynamics of open system in the pure dephasing models 
with the boson and spin baths. We found that the ini- 
tial correlation can reveal the memory of bath during the 
time evolution. The decoherence of a qubit coupled to a 
boson bath is more enhanced with respect to a spin bath 
in the short time. We also demonstrated that the trace 
distance between two states of qubit can increase above 
its initial value witnessing the initial non-classical corre- 
lation between the qubit and bath, and that the entan- 
glement of two qubits, locally interacting with an inde- 
pendent bath, can damply oscillate and revive to a large 
amount compared to the uncorrelated initial condition. 
We finally showed that the initial correlation affects the 
crossover of decoherence from the dynamical enhance- 
ment to suppression under the bang-bang control pulses 
to some extent only at the beginning of evolution, and 
induces more changes for the boson (weak) bath than for 
the spin (strong) bath. On the other hand, the higher 
temperature for the boson bath, the more efficient for 
the control pulses to make the crossover occur. For the 
spin bath, the crossover under control shows the temper- 
ature insensitivity instead. 
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